Shards of Broken Symmetry: Topological Defects as Traces of the Phase Transition 
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We discuss the origin of topological defects in phase transitions and analyze their role as a "di- 
agnostic tool" in the study of the non-equilibrium dynamics of symmetry breaking. Homogeneous 
second order phase transitions are the focus of our attention, but the same paradigm is applied to 
the cross-over and inhomogeneous transitions. The discrepancy between the experimental results 
in 3 He and 4 He is discussed in the light of recent numerical studies. The possible role of the 
Ginzburg regime in determining the vortex line density for the case of a quench in 4 He is raised 
and tentatively dismissed. The difference in the anticipated origin of the dominant signal in the 
two ( 3 He and 4 He) cases is pointed out and the resulting consequences for the subsequent decay 
of vorticity are noted. The possibility of a significant discrepancy between the effective field theory 
and (quantum) kinetic theory descriptions of the order parameter is briefly touched upon, using 
atomic Bose-Einstein condensates as an example. 



INTRODUCTION 

The theory of the creation of topological defects ap- 
peals to models of critical dynamics and to our under- 
standing of the processes which occur when phase tran- 
sitions take place. Consequently, topological defects can 
be used as "symptoms", macroscopic manifestations of 
underlying physical processes, which in turn can help di- 
agnose the nature of critical dynamics. 

For first order phase transitions there is little doubt 
that nucleation - a process understood for over half a 
century - is an essentially accurate, universal yet sim- 
ple model. A similarly simple model of the dynamics of 
second order phase transitions was proposed much more 
recently [H-0|. One of its implications is the ability to 
predict the size of the ordered patches of the new lower 
symmetry phase, right from its inception. This allows 
one to calculate the initial density of topological defects 
through the estimate put forward in the seminal paper 
by Kibble 0, HI . It may also lead to a revision of the sce- 
narios for baryogenesis and chiral symmetry restoration 

0] ) as well as other related phenomena (such as the 
A-B 3 He transition, see [1, Some of the predictions 
based on the new paradigm have been successfully tested 
and refined in numerical experiments [Ml. 

More im- 
portantly the prediction of copious vortex production in 
superfluid phase transitions has been experimentally ver- 
ified in 3 He by two very different strategies in two distinct 
parameter regimes (l5l.[l6|. The situation in 4 He [l7| and 
the initial indications from high temperature supercon- 
ductors [l8j are, however, at best inconclusive. Indeed 
there are still differences concerning analytic estimates 
of the initial density of defects in the underdamped case 
[l9l [2fj| . even in ID. The aim of this paper is two-fold: 
We shall start with a brief summary of the paradigm on 
which the emerging understanding of second order phase 
transitions is based. We shall then explore its extrapo- 



lations and investigate the experimental, numerical and 
analytic evidence for and against this mechanism, in var- 
ious settings. This paper is not really an introductory 
survey to the extent to which our lectures were. We have 
decided that the existing literature (including the reviews 
of Zurek [5H and Eltsov, Krusius and Volovik [52| as well 
as the other papers mentioned above and below) already 
serves this purpose. Rather we aim to perform a "re- 
connaissance by force" of what is likely to be the most 
interesting "proving grounds" for the ideas summarized 
briefly in the following section. 



CRITICAL DYNAMICS AND DEFECT 
FORMATION 



Second order transitions fall into universality classes 
which are characterized by the behavior of the healing 
length £ and the relaxation time r (among other quanti- 
ties) as a function of the relative temperature 



e(T) 



T-T c 



(1) 



Thus, t ~ \e\~ uz and £ ~ | e | ^ diverge in the vicin- 
ity of e = 0, where v and 2 are universal critical expo- 
nents. A very specific model which represents a large 
class of second order phase transitions is the so-called 
Landau-Ginzburg theory There, the dynamics of the or- 
der parameter is thought to effectively obey a Langevin 
equation of the form: 

(p + rjip - c 2 S7 2 ip + i [Xip 3 + m 2 4>] = 0(t, x). (2) 

Above, r\ characterizes the viscosity in the system, 
while c, A are constant coefficients. The mass term m 2 
can depend explicitly on time. Moreover, the correlation 
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function for the noise 

(0(t, x)0{t', x')) = 2f]QS{x - x')5(t - t'), (3) 

includes the temperature parameter 0, which can vary. 
The change of m 2 , eg. m 2 = mge(i) or of or both, may 
precipitate the phase transition. 

We shall assume that in the vicinity of the critical tem- 
perature e(t) obeys a simple relation 



e(t)=t/r Q . 



(4) 



In that case the dynamics of the order parameter can 
be approximately divided into the adiabatic and impulse 
regimes . with the boundary which occurs at time t 
when the relaxation time of the order parameter equals 
the characteristic time on which e(i) changes: 



r(e(i)) 



e(t) 



(5) 



The timescale on which the order parameter will react 
to changes of e(t) depends on whether <p or (p dominates. 
For the Landau-Ginzburg theory, in the two cases 



| e |2 : 



(6) 



respectively, where tq — 1/rao, and mo is the mass term 
in Eq. ©, evaluated for T = (i.e. for e = -1). 

Using e(i) = £/tq, we can now solve for t, to obtain: 



tu 



, 2/3 1/3 



(7) 



To estimate the scale of the domains which could have be- 
come uniform through dynamics in the adiabatic regime, 
we should need: 




= ± 



2/3 



(8) 



The characteristic scale is then given by £ = £o / 
which yields: 



1/4 



£<p — £0 — 2 



TO 



1/3 



(9) 



The initial density of defects can now be estimated using 
an argument due to Kibble Q, which will imply a £-sized 
unit of defect per £ sized volume, i.e. n ~ l/£ for vortex 
strings in two spatial dimensions. 

The above calculation is based on the assumption that 
the order parameter approximates its equilibrium con- 
figuration until — t, at which point it ceases to evolve 
dynamically (although noise and damping continue un- 
abated). The dynamical evolution restarts at +t, below 
the critical point, but by then it may be too late to undo 



non-trivial topological arrangements of <p inherited from 
above T c . 

This same paradigm decides when the overdamped or 
underdamped estimates are relevant. For, in view of the 
above argument, it is essential to decide whether the dy- 
namics of the order parameter is overdamped at t, i.e. 
whether; 



(10) 



This can be evaluated directly from Eq. © with the help 
of the above estimates for t and e, and leads to the in- 
equality: 



(777-0) >t /t q . 



(11) 



Numerical studies have by now confirmed this 
paradigm. The scalings which we obtained follow the- 
oretical predictions both when the quench is induced by 
an explicit change of the mass term in Eq. ([2]) (lol 11 1 , (l2l | 
and when the temperature of the noise is adjusted, but 
m set to a constant [13j. Moreover, the switch from the 
overdamped to the underdamped behavior occurs where 
expected, and with the consequences consistent with the 
scaling implied by the paradigm Eqs. (|8|)- ([TT1) [IOL lll|. 
fl2l | . The same reasoning can be of course repeated using 
other values of critical exponents relevant for other cases 
f23l | , which has been already done in some cases [2l[ , [TJ[ • 

While the scalings accord well with the theoretical pre- 
dictions, the specific density of defects n is lower than the 
appropriate inverse power of £ ; 



1 



(12) 



where / is always more than unity, and usually in the 
range 8-15 [Tj-dJ. We shall return to its estimates 
later in this paper. 



CROSSOVER TRANSITIONS 

An interesting case of transitions which does not con- 
form to the second-order universality class occurs when 
the critical scaling behavior of the healing length and 
of the relaxation time "tapers off" (i.e. is fully analytic) 
very near to e = 0. For instance a hypothetical relaxation 
time and healing length dependences; 



To 



A)' 



£0 



(13) 



illustrate such a crossover transition. 

Examples of crossover phenomena are ubiquitous. One 
of the most interesting cases comes from the study of the 
electroweak standard model, where, for Higgs masses not 
yet excluded by experiment, the transition appears to be 
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a crossover [24|. A crossover transition may also occur 
in the presence of impurities, anisotropics, weak external 
fields, or finite size scaling, instead of the expected critical 
behavior of the ideal model. It also substitutes critical 
behavior when non-perturbative fluctuations exist in the 
spectrum of the theory which are favored entropically 
and can destroy long range order. An example of the 
latter is a \<j) A theory or a (short-range) Ising model in 
one spatial dimension. 

Thus, t and £ in the one-dimensional cases investi- 
gated numerically (lol . 11 1 as examples of the second order 
phase transition are expected to taper off in the imme- 
diate vicinity of the critical temperature. Presumably 
this occurs for very small values of A and S, so that the 
scaling behavior encountered in the vicinity of e is not af- 
fected. Nevertheless, it is interesting to investigate what 
does the paradigm predict in the case of such crossover 
transformations. 

We follow the footsteps of the argument outlined in 
the preceding section, and obtain t by solving : 

r(e(t)) = t, (14) 

which now leads to the quadratic equation 

\t\ 2 + \1\Atq - tqTq = 0. (15) 

Consequently 



- At q + t q v/A 2 + 4r /r Q 



and 



k"l = -y- + ^A 2 +4ro/r Q , 



(16) 



(17) 



where we have picked the physically relevant root of 
Eq. (USD. 

We note that in the limit of a 'real' second order phase 
transition (A — > 0) we recover the old result, Eq. (J5J, 
providing that the change of notation (to now used to be 
TyTg in Eq. (8])) is acknowledged. On the other hand, when 
the quench is very slow and A 2 >> 4to/tq, e — > Ato/tq, 
which itself is small compared with A (and, presumably 
also S since A ~ 5 can be expected). Consequently for 
relatively rapid quenches 



(18) 



which approaches Eq. © for sufficiently small S. In the 
other limit; 



(19) 



and the size of the coherent domains of the order param- 
eter saturates. 



We note that the above discussion should be regarded 
more as an exercise in extending the paradigm rather 
than as a generically valid theory, applicable to all 
crossover phase transitions. In particular, in some cases 
second order transitions may change into crossovers when 
an external bias which influences the choices of the bro- 
ken symmetry vacuum is introduced. In such cases the 
externally imposed (rather than spontaneous) symmetry 
breaking will favor a particular vacuum and will lead to 
a suppression of topological defect production [2(|, Q). 

Moreover, in case of the crossover transitions the in- 
fluence of the Ginzburg regime may need to be carefully 
examined as its role in the generation and survival of 
topological defects is still a subject of dispute. 



INHOMOGENEOUS TRANSITIONS 

Homogeneous quenches are a convenient idealization 
and may be a good approximation in some cases. How- 
ever, in reality, the change of thermodynamic parameters 
is unlikely to be ideally uniform: 

1) Experiments carried out in 3 He [Tilllij. where a small 
volume of superfluid is re-heated to normal state, and 
subsequently rapidly cools to the temperature of the sur- 
rounding superfluid, are a good example of an inhomoge- 
neous quench: The normal region shrinks from the out- 
side. Yet, topological defects are created, thus suggesting 
that the phases of distinct domains within the re-heated 
region are selected independently. 

2) Another example are relativistic heavy ion collisions 
where, according to Bjorken scenario [25], a finite vol- 
ume of quark-gluon plasma can be created. The plasma 
expands in the direction of collision and cools from the 
outside in the perpendicular direction. The phase transi- 
tion in this case can be first or second order (or a smooth 
crossover) depending on the parameters of the collision. 

3) Any generic experiment based on pressure and/or tem- 
perature quench is to some degree inhomogeneous be- 
cause of finite velocity of sound and finite heat conduc- 
tance. 

The mass parameter e(t, r), varying in both time and 
space, must be considered in defect formation. As a con- 
sequence, locations entering the broken symmetry phase 
first could communicate their choice of the new vacuum 
as the phase ordered region spreads in the wake of the 
phase transition front. When this process dominates, 
symmetry breaking in various, even distant, locations is 
no longer causally independent. The domain where the 
phase transition occurred first may impose its choice on 
the rest of the volume, thus suppressing or even halting 
production of topological defects. This happens if veloc- 
ity of the critical front is less than certain characteristic 
velocity. 
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Second Order Transition 



where 



The characteristic velocity in an overdamped transi- 
tion can be estimated as follows: The freeze-out healing 
length is set at t as £ = £o (tq/to) 1 ■ At the same in- 
stant the relaxation time is f = (tqTq) 1 / 2 . These two 
scales can be combined [2j to give a velocity scale 



(20) 



where v = W T o- 

The density of defects ./V as a function of critical front 
velocity is expected to change qualitatively at v. Above 
v the homogeneous estimates should hold. Below -0 the 
density should be suppressed. Kibble and Volovik [26] 
suggested that TV ~ v/v for small v < v. Dziarmaga, 
Laguna and Zurek [27| argued that N is exponentially 
suppressed below v. There is qualitative difference be- 
tween the two proposals. The former option suggests 
that however one makes a quench one will always get 
some defects, the latter implies that if one's inhomoge- 
neous quench is sufficiently slow one will get no defects 
at all. In what follows we will quantify what "sufficiently 
slow" means. 

Decay of the False Vacuum As a simple warm up 
exercise, let us consider decay of a false symmetric vac- 
uum to a true symmetry broken ground state in a one- 
dimensional dissipative ip 4 model 



d$ (f + r]d t tp 



e<p) = , 



(21) 



where ip(t, x) is a real order parameter and e measures 
the degree of symmetry breaking i.e. m 2 = — e. Without 
loosing generality, we look for a solution ip(t, x) which 
interpolates between tp(t, — oo) = — y/e and tp(t, +oo) = 
0. Such a solution can not be static. It is a stationary 
half-kink 



<p(t,x) 



-\/e 1 + cxp 



\/e(x- v t t) 



2 y/T=V. 
moving with characteristic velocity 



-1/2 



jj-s-cx) 3\/e 
~2V 



(22) 



(23) 



It is worth noting that the decay velocity vt increases 
with e. 

Shock Wave Our shock wave inhomogeneous quench 
model consists of a sharp "pressure front" propagating 
with velocity v; that is, 



%<p+ndttp- d 2 x <p+-(<p 3 -e(t,x)cp) = 0(t,x) , (24) 



e{t,x) = Sign(£ — x/v) 



(25) 



is the relative temperature and 0(t,x) is a Gaussian 
white noise of temperature with correlations given by 
Eq. ©. 

There are two qualitatively different regimes: 

1) v > Vt, the phase front propagates faster than the 
false vacuum can decay. The half-kink (f2"2")l lags behind 
the front (1251) ; a supercooled symmetric phase grows with 
velocity v — v t . The supercooled phase cannot last for 
long; it is unstable, and the noise makes it decay into the 
true vacuum. 

2) v < v t} the phase front is slow enough for a half- 
kink to move in step with the front, tp(t,x) — H v (x — vt) . 
The symmetric vacuum decays into one definite non- 
symmetric vacua. The choice is determined by the 
boundary condition at x — > — oo. No topological de- 
fects are produced in this regime. The stationary solution 
H v (x — vt) is stable against small perturbations [2?J ■ 

These expectations are borne out by the numerical 
study of kink formation in [27]. Numerical results are 
presented in Fig. 1. 

Linear Front Let us consider now a system in which 
the inhomogeneous quench takes place via linear transi- 
tion 



e(t,x) = (t-x/v)/T Q 



(26) 



In the absence of noise, the propagating front is fol- 
lowed by a stationary half-kink. This half-kink moves 
somewhat behind the front. Its location is determined 
by the place where the threshold velocity (|2"3")) is equal to 
the front velocity, v t [e(t,x)] = v. The distance between 
the front and the half-kink increases as v 3 . This distance 
gives the size of the supercooled region. When the su- 
percooled region is narrow then it is stable against small 
perturbations so that no defects are produced. If 



v > v t 



11.7 



-1/2 



T£» 3^ ( JL)i/4 s 4.07 v 

V TQ 



(27) 



then the region is broad enough to be unstable [27| and 
the production of defects is no longer suppressed. 

This prediction is confirmed by the numerical study of 
linear quenches in Ref. |27j, compare Fig. 1. However, the 
threshold velocity apparently gradually decreases with 
increasing noise temperature 0. This decrease of the 
threshold for kink formation is due to the thermal nucle- 
ation of kinks. Quantitative estimates for this effect are 
given in (2f| . 
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Figure 1: a. Density of kinks n as a function of velocity v 
for the shock wave 1|25[) with rj = 1 (overdamped system). 
In this overdamped regime, the predicted threshold velocity 
is vt = 0.83. The plots from top to bottom correspond to 
9 = 1CT 1 , 10" 2 , 10" 4 , 10" 6 , 10 -8 , 10 -10 . At low 8, we get 
a clear cut-off velocity at v « 0.8, which is consistent with the 
prediction, b. Density of kinks n as a function of velocity v for 
the linear inhomogeneous quench, Eq. (I2T>|) , with tq — 64 and 
rj = 1. The predicted threshold is v t = 0.77. This cut-off is 
achieved for low Q. The plots from top to bottom correspond 
to 9 = 10 _1 , 10~ 2 , 10" 4 , 10~ 6 , 10~ 8 , 10~ 10 . 



First Order Transition 

We assume the transition is strongly first order and 
that it goes by bubble nucleation. To be more specific 
we consider a toy model in 3 dimensions 



where tp. is real order parameter. The effective potential is 
of the tp 6 type. Provided that b 2 > 4ac, it has symmetric 
minimum at tp — and two symmetry broken minima 

at tp = ±ip m = ±\J (b + \Jb 2 — Aac)/2c. We assume that 
b, c are constant and that symmetry breaking transition 
is driven by a decreasing below its critical value a c = 
3b 2 /16c. At a = a c all three minima are degenerate. 

Decay of the False Vacuum Suppose that a < a c . 
Let us consider decay of the false symmetric vacuum to 
the true symmetry broken phase in a one dimensional 
version of the model Eq. (|2"51) . We look for a solution 
which interpolates between tp — tp m for x — > — oo and 
ip = for x — > +oo. The solution is a stationary half- 
kink H(x — vtt) moving with velocity 



-b + 2y/b 2 - 4ac 



which has an envelope function 



H(x) 



^7) 



CXp Q;f 

2c 



(29) 



(30) 



where a — y^Ac/3ip 2 n . This way the false ip = vacuum 
decays into the true tp = tp m vacuum in the absence of 
noise. The decay velocity vt is zero for a = a c , it increases 
with increasing supercooling or with decreasing a. 

Shock Wave In the shock wave model a sharp front 
propagates with velocity v 



a = a c — Aa Sign(i — x/v) . 



(31) 



Similarly as for second order transitions there are two 
regimes: 

1) v > Vt, the pressure front propagates faster than 
the false vacuum can decay. The half-kink lags behind 
the front. The supercooled phase in between them grows 
linearly with time. The phase is unstable, it decays by 
bubble nucleation just as for a homogeneous transition. 
Homogeneous estimates of defect density apply in this 
case. 

2) v < vt, the half- kink is faster. It moves in step 
with the front while its tail penetrating into the symmet- 
ric phase. There is no supercooled phase where bubbles 
could be nucleated. The symmetric phase goes smoothly 
into one of the symmetry broken phases. 

Linear Front Let the inhomogeneous quench pro- 
ceed by a linear front moving with velocity v 



= a c -(t- x/v)/t q . 



(32) 



d t tp = \7 2 tp — atp + btp 3 - ctp 5 



O 



(28) 



The half-kink follows the critical front staying at a cer- 
tain distance behind it. The distance D is such that the 
half- kink velocity vt, which depends on the local value 
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of a, is equal to the front velocity v, Vt(a) — v. With 
increasing v the half-kink settles at increasing values of 
local a. Close to the critical front the radius of the criti- 
cal bubble is infinite and at the same time the nucleation 
rate is infinitely small. As we go away from the front in 
the direction of the half-kink the critical radius shrinks. 
At a certain distance L from the front the energy of the 
critical bubble becomes comparable to the temperature 
O. At this point bubble nucleation becomes possible. If 
L < D bubbles can be nucleated in the supercooled re- 
gion between the front and the half-kink. If L > D then 
there is no bubble nucleation and no defects can be born 
in the supercooled area. 

Now we estimate the critical velocity such that L = D. 
The half-kink is located at such an a that v t (a) = v. 
L = D providing that for this a the energy of the critical 
bubble E{a) is equal to temperature 6. The critical bub- 
ble is a metastable spherically symmetric static solution 
of Eq. (|28p with, say, ip m vacuum inside and vacuum 
outside its wall. Its energy can be easily estimated when 
the width of its wall is negligible as compared to its radius 
R c (a). An approximate solution is given by H[r — R c (a)], 
where the critical radius is 

R c (a) = V ^ . (33) 

V ; -b + 2Vb 2 - Aac 

The energy of the critical bubble E{a) has a nega- 
tive volume contribution, (4irR3/3)V(tp m ), and a posi- 
tive surface tension term, (AirR 2 ) J dx \H'(x)] 2 . When 
the solution of v t (a) = v is put into E(a) and then the 
equation E(a) — is solved, one obtains a critical veloc- 
ity 

V "={ 4 C 36 -) (34) 

for L = D. For v > v cr bubbles can nucleate in be- 
tween the half-kink and the front and thus the necessary 
condition for topological defects production is satisfied. 

The formula for v cr , Eq. (|34p . is still a crude lower es- 
timate for the critical velocity. In fact it is not sufficient 
to nucleate some bubbles. Individual bubbles would co- 
alesce with the half-kink without any chance to trap any 
nontrivial winding number. The bubbles should be nu- 
cleated in large numbers or have enough time to grow so 
that they can mutually coalesce before merging with the 
half-kink. Still, the argument which leads to v cr demon- 
strates that there is a threshold velocity for defect for- 
mation. 

Higher Dimensions 

The theory can be generalized to higher dimensions 
and to a complex order parameter in a straightforward 



manner. Its major result is that a subthreshold inhomo- 
geneous quench does not produce any variation of the or- 
der parameter in the direction normal to the front. This 
excludes any possibility of production of vortex loops or 
closed membranes entirely contained in the bulk, as well 
as of any pointlike defects. Some extended defects can 
grow into the bulk provided their seeds were created at 
this edge of the system where the symmetry was bro- 
ken first. In first approximation such, say, vortices grow 
into the bulk, following the passing front, while keep- 
ing their direction normal to the front. In the end we 
do not get any chaotic tangle of strings and string loops 
but parallel "combed" vortices. There are two important 
perturbations to this "combed" picture: 

1) Thermal fluctuations make the strings look more 
random but without backtracking and with string ten- 
sion tending to smooth the small scale fluctuations. The 
ends of the strings and antistrings at the critical front 
are wandering around. Eventually an end of a string and 
of an antistring may meet so that the strings join into a 
half-loop with its both ends attached to the initial edge 
of the system. String tension shrinks the half-loop to the 
edge where it unwinds. 

2) A much more efficient factor to remove vortices from 
the bulk are their mutual interactions. Global parallel 
string and antistring attract one another so that their 
ends at the critical front do not seek each other at ran- 
dom but tend to fuse in a deterministic way. This mecha- 
nism makes the number of strings in the bulk decay with 
increasing distance between the front and the initial edge. 

The factors (1) and (2) lead to a picture in which the 
critical front initially draws some parallel strings and an- 
tistrings from the edge, then the strings recombine by 
joining ends and shrinking back to the edge. In the end 
only the net surplus of strings (or antistrings) is left in 
the bulk. 

These ideas are supported by experiments: 

1) Disclinations produced during a quench from disor- 
dered to nematic phase in liquid crystals [29]. This is a 
weakly first order transition. In early attempts to make 
cosmological experiments in liquid crystals the disclina- 
tions were observed to grow approximately combed, join 
ends and shrink to the initial edge. Later on it was real- 
ized that these quenches were not homogeneous enough 

2) Czochralski method of growing monocrystals, which 
is widely used to grow silicon monocrystals necessary for 
microchips. In this method, discovered in the thirties, a 
surface of liquid material is touched with a monocrystal 
template. As the template is slowly lifted up it drags 
a column of crystal out of the container. The top part 
of the column is cold while its bottom part is at the 
melting temperature - the transition is inhomogeneous. 
If the template is lifted slowly enough, then no defects 
of the crystal lattice are produced which might spoil the 
monocrystal. 
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To conclude this section: in an inhomogeneous quench 
there is a threshold velocity Vt of the critical front. Above 
the threshold defects are produced like in a homogeneous 
quench. Below the threshold one gets no defects; instead 
a clean monocrystal or a "disoriented chiral condensate" 
is grown with a vacuum which may be uniform over signif- 
icant distances, but which differs from the true vacuum. 



DEFECT FORMATION AND THE GINZBURG 
REGIME 

Recently a new 4 He experiment [13] was devised, im- 
proving on the apparatus used earlier by McClintock et 
al. [31] to implement a superfluid transition in 4 He 
through a sudden pressure quench. The corresponding 
results are rather surprising. They show no evidence for 
the formation of topological defects at the anticipated 
levels, contrary to expectations based both on the old ex- 

3Hl- 



periment [31], the theory |42| and the 3 He data 
The discrepancy with the earlier 4 He quench data is now 
seen as the evidence of mechanical stirring in the first 
version of the experiment. Nevertheless to address the 
discrepancy with 3 He it was suggested (32j that because 
the Ginzburg regime in 4 He extends over a broad range of 
temperatures around the A-line, large scale fluctuations 
may be able to unwind and alter the configuration of the 
order parameter (in contrast to 3 He) while the quench 
proceeds. 

The Ginzburg temperature is defined through the loss 
of ability of the order parameter to hop, through ther- 
mal activation, over the potential barrier between broken 
symmetry vacua. Thus one might worry with Karra and 
Rivers [32| that when the defect densities are eventually 
measured, at a much later time, little or no string would 
have survived unwinding through thermal activation. In 
this section we investigate this possibility and more gen- 
erally report a numerical study of the effect of thermal 
fluctuations on topological defect formation and evolu- 
tion. 

Originally the Ginzburg temperature T G was suggested 
to be the time of formation of topological defects 
since, at lower temperatures, thermal fluctuations would 
be unable to overcome the potential energy barrier asso- 
ciated with the defect's topological stability. 

In reality the situation is more complex. In equilib- 
rium at any given temperature T (including of course 
temperatures in the Ginzburg regime) a range of string 
configurations will exist. However, long strings can only 



33] 



exist in equilibrium strictly above T c 

To freeze them out, i.e., to form them, energy (asso- 
ciated with the string tension) must be extracted from 
the system. This necessarily breaks time invariance and 
will lead to locally preferred nonequilibrium field config- 
urations. Subsequently the system will order over larger 
and larger spatial scales, leading to mutual string anni- 



hilation. 

The initial density of defects entering this stage of evo- 
lution is computed by the theory of section II. This den- 
sity is set at an effective temperature — e, which in A He is 
well within the estimates for the width of the Ginzburg 
regime. By contrast, in 3 He the Ginzburg temperature 
is small compared to the typical e. What happens to 
the initial densities of string when the system is exposed 
to temperatures in the Ginzburg regime for an extended 
amount of time? 

In order to investigate this issue we need a quantita- 
tive definition of Tq. In tune with the arguments given 
above consider a volume of characteristic size £(T), the 
correlation length, and a theory with two energetically 
degenerate minima of an effective potential V(<p), sepa- 
rated by a potential barrier AV. The rate for the field 
to change coherently from one minimum to the other per 
unit volume due to thermal activation is exp [— AV/fesT]. 
For an effective potential of the form (obtained, eg. per- 
turbatively at 1-loop) 

V(0) = -im 3 (ry + A V 4 (35) 

AV = TO ^p . For a volume £ 3 , we define T G such that 
the probability of overcoming the potential barrier is of 
order unity: 

To : ^p^ 3 (T G ) = 1 «• XT G /m(T G ) = \ . (36) 

This definition however has some caveats, for instance, 
an effective potential of the form Eq. ([33)) is only valid 
for the mean field and not on smaller scales. A more 
careful accounting of scales leads to different results [lH , 
which show an enhancement of the hoping probability. 
Thus, the factor of 1/4 in Eq. (f3"6"| should not be taken 
at face value. 

A more rigorous definition arises from the range of tem- 
peratures below T c for which fluctuations are large and 
consequently where perturbative finite temperature field 
theory fails to be useful. In order to set up a pertur- 
bative scheme at finite temperature from an initial 3+1 
dimensional quantum field theory one implements dimen- 
sional reduction which is valid provided the temperature 
is high compared to all mass scales. As a consequence 
the coupling of the dimensionally reduced 3D field the- 
ory becomes dimensionful, i.e. A — > AT = A3. In order 
to proceed one has to identify an appropriate dimension- 
less coupling. This is done by taking XT/m(T). The 
Ginzburg regime is entered when this 3D effective cou- 
pling becomes strong, in the vicinity of the critical point, 
namely 



\T G /m(T G ) = 1. 



(37) 



To compute T G one needs the scaling of m(T) in the 
critical domain. We write m 2 (T) = rn^e", with e being 
the reduced temperature e 
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Figure 2: Temperature trajectories for testing the effect of ex- 
posure to the Ginzburg regime on string densities. The system 
is first thermalized at a high temperature and then placed in 
contact with a heat bath at an intermediate temperature e; 
below T c , for a time interval At. 



Figure 3: The string density measured at a later time t » At 
vs. intermediate temperature e*. From top to bottom the 
three plots correspond to At = 10, 20, 50, during which the 
system remained in contact with a heat bath at Ti. There 
is no visible role played by intermediate temperatures within 
the Ginzburg regime. 



Thus 6q — —0.18 for v = 0.5. This mean-field estimate 
produces an upper bound in T for Tq (and lower bound 
for j3 — 1/T). For realistic 3D exponents, v — 0.67, we 
obtain ec = —0.25. The first criterion, based on the hop- 
ping of a correlation sized volume, results in higher values 
of Tq. This brings about a relatively large uncertainty 
in the value of Tq, which is 18 — 25% below T c . 

Strings Survive the Ginzburg Regime 

In order to investigate the role of the Ginzburg tem- 
perature in the dynamics of defect formation we deliber- 
ately expose the system to a heat bath at temperature 
ej, within the Ginzburg regime and below. We repeat 
this procedure for a range of time intervals At, after 
which the bath temperature is taken to zero. This set 
of temperature trajectories is shown in Fig. 2. We are 
attempting to emulate the worst case scenario of an ex- 
perimental quench where the temperature or pressure are 
dropped monotonically but where the system makes a 
long stopover within the Ginzburg regime. 

We would expect that, if the Ginzburg regime indeed 
produced enhanced decay of strings, then the string den- 
sities measured at later times should be smaller the longer 
the time the system spent within the range T c > T > Tg- 

We have measured the final string densities at a time 
t » At. 

Our results for the final string densities as a function 
of intermediate temperature ej and At are summarized 



in Figure 3. There is no apparent effect of the Ginzburg 
regime in reducing string densities at formation. 

If any trend is visible from Figure 3 it is the opposite, 
namely that the lower e», the less string is measured at 
later times. This is consistent with the relaxation of the 
string network, resulting in vortex annihilation controlled 
by the string tension (which is smaller near T r ) and with 
results the thermodynamics of vortex strings [33] . 

Memory of the Order Parameter Configuration near 

T c 

An independent test on the possible role of thermal 
fluctuations in affecting string densities consists in re- 
heating a quenched system to a temperature around T c 
(both below and above it) and cooling it again. This 
process tests the memory of the order parameter as well 
as that of other related quantities [see also [Hj]], such as 
defects. These temperature (e(t)) trajectories are illus- 
trated in Fig. 4a. 

We are particularly interested in investigating under 
what circumstances thermal fluctuations can affect the 
large scale configuration of the order parameter. 

We define the unequal-time correlation function 

(ip{x,t lh )ip(x,t + f rh )) = 

2 

N/l/lVj i x i i tih)<Pj {xj,t + t rh ), (38) 

j=l i 
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Figure 4: a. Dependence of the bath temperature e in time. 
After being quenched in temperature (tq = 80) the sys- 
tem is reheated at the same rate to a temperature e/ = 
0.469, 0.256, 0.061, -0.068 (top to bottom) 'and cooled again, 
b. The correlation function between the field at the time just 
before reheating and at later times, (y>i(i r h, x)tfii(t + t^jx)) 
is plotted. There is a universal short time transient for the 
decorrelation of the field over small scales while the long time 
tails of the correlation function describe change over the mean 
fields. All four trajectories cross the Ginzburg regime, but 
only those reaching or crossing +e display a significant mem- 
ory loss. 



where N is an irrelevant normalization factor. Note 
that (p(x, t r h)(p{x, t + i r h)) cannot be complex as we have 
summed over the field's components. This correlator has 
several interesting properties. For short times it displays 
a characteristic time, which describes the decay of corre- 
lations over very small spatial scales. This is the initial 
transient in Fig. 4b. For later times the residual correla- 
tion comes from the motion of the order parameter (the 
field volume average). This average can be either posi- 
tive or negative but, if thermal, will converge to zero at 
and above T c . 

Now, we are interested in determining whether the final 
field configuration over large spatial scales is correlated 
to the configuration prior to reheating. Fig. 4 shows that 
only if one crosses T c , by more than +e, is the memory 
of the initial quenched configuration erased (see in par- 
ticular the two trajectories reaching higher temperatures 
in comparison to the others). For these trajectories the 
field correlations reach zero and after reheating evolve to 
a value manifestly different from that prior to reheating. 

For trajectories within the Ginzburg regime, that do 
not cross T c , the change in the configuration of the order 
parameter as measured by Eq. (|38l) is small. In partic- 



ular the field configuration existing before reheating is 
approximately recovered as the fields are cooled. The 
same is true for the string densities, including those of 
long strings. 

Thus we are led to conclude (see also [33|) that even 
prolonged exposure of a quenched field configuration to 
the Ginzburg regime has little consequences in changing 
the order parameter configurations emerging at — e, and 
associated string densities. In addition we have shown 
that to truly destroy a quenched field configuration ex- 
isting below — e, one has to expose the system to tem- 
peratures well above T c . In particular for any particular 
quench trajectory, a temperature of T ~ T c + e, must be 
reached and maintained for a time ~ t in order to erase 
memory of the initial configuration. 

These results fully support the theory of section II for 
the critical dynamics of second order transitions and all 
known thermodynamic results for vortex strings in O(N) 
theories. Thus we expect the results of this section to 
carry over from our models to the Lancaster He exper- 
iments. The results of Ref. [l7j in these experiments 
cannot therefore be attributed to the effects of Ginzburg 
regime in 4 He. 

In the next section we offer an alternative explanation. 



WHAT IS BEING OBSERVED IN THE 3 HE AND 
4 HE EXPERIMENTS. 



The several experiments in Helium, and more recently 
in superconductors, testing the theory of defect forma- 
tion rely on substantially different processes to induce 
the phase transition and measure defects. 

In this section we analyze, in the light of our own the- 
oretical results, how experimental procedures can lead to 
the detection of substantially different defect densities. 

Two particular factors play a decisive role in the value 
of the topological defect density measured - the time and 
procedure of measurement after the quench and the ini- 
tial/final state of the system. 



The Lancaster experiments in He 

In the Lancaster experiments in 4 iJe the defect density 
is measured through the attenuation of a second sound 
signal (a heat pulse). This probe can only detect den- 
sities above a certain threshold (if the theory of section 
II is used / < 10 would be required, which is at odds 
with the results of the numerical studies (lol . 12 1 and es- 



pecially [13j|). Moreover, the density at formation has to 
be extrapolated from the data obtained at relatively late 
times - the signal is noisy shortly after the quench [l?} • 

After being formed by the critical dynamics of the 
phase transition vortex strings decay away, as the sys- 
tem orders and cools. This decay has been modeled by 
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+t and the "time of formation" ({\<j>\) = 0.95), for tq = 64. 
Data sets denoted by (*,A,D, •) correspond to increasingly 
later times. Lines show the integral distributions, eg. L.(l) = 

Figure 5: The decay of vorticity under a quench. Initially the ^ ( °° I' Ld.(l')dl' . It is clear that at late times only long strings 

defect network includes both long strings and small loops. At survive. 

late times the network coarsens and only long strings remain, 

see also Fig. 6. 



Vinen's Equation 

n=- 7 n 2 , (39) 

where n is the length of string per unit volume, i.e. the 
string length density, 7 > 0. This model has been ob- 
served in the same experiment to describe very well the 
decay of vorticity induced initially through a fluid flow. 

Vorticity created thermally is potentially different from 
that formed under an external flow. We know from sev- 
eral theoretical and numerical indications that a thermal 
distribution of vortices close to the transition is com- 
prised of both long strings and small loops, see Fig. 5. 

These two populations decay very differently in the 
wake of the quench. Without any mechanism for stabi- 
lization the loops tend to disappear in a fast transient. In 
contrast the long strings loose some of their small scale 
structure but survive, and will ultimately set the decay 
pattern described by the Vinen equation. 

It is the surviving long strings - eventually measured - 
that will provide the experimental signal in the 4 ffe ex- 
periments. This is shown in Figs. 5 and 6. As the system 
is quenched from higher temperatures or pressures, an 
initial string network comprised of strings of all lengths 
looses its loops and settles to a much slower decay trend 
dominated by long strings. 

The crucial question then is whether enough long 
string would persist at the time of measurement to yield 



a positive signal. 

We have performed a very similar procedure in a nu- 
merical "experiment" [l3j]. We observed that at long 
times string densities could be measured that agreed ex- 
tremely well with all features of the theory. Our defini- 
tion of long times was intimately connected to the com- 
pletion of the phase transition expressed in the expecta- 
tion values of the order parameter (|<^|) ~ 0.9 — 0.975. 
The effective / measured then was in the range / 
II Ifi. All our indications are that the Lancaster He 
experiment performs its measurements much later (up to 
2 orders of magnitude) than we do, thus leading to even 
smaller string densities. Such string densities could evade 
detection under the second sound ex per imental probe, 
which may lead to the negative result |17| . 

The other possibility is that the annihilation of a net- 
work of vortices arising from a rapid quench may proceed 
at a rate different than for the network produced by tur- 
bulent flow. Clearly, in case of turbulent flow there may 
be correlation between the orientation of nearby vortices. 
By contrast, vortices created by quench are anticorre- 
lated (see, e.g., discussion in section 4 of Ref.(2l]|). It 
seems plausible that annihilation of vortices would be 
faster in that case. In order to measure a positive signal 
in these circumstances the measurements would have to 
be made sooner, after a much faster quench, or with a 
higher sensitivity. 
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Experiments in 3 He 

In contrast to the experiments in He described above 
which appear sensitive to the "infinite" string, experi- 
ments in 3 He have pursued two independent strategies 
both of which allow one to stabilize and detect defect 
loops of various sizes: either vorticity is stabilized and 
amplified by the flow and then measured directly using 
nuclear magnetic resonance (the Helsinki experiment) or 
it is inferred from a certain amount of missing energy 
(the Grenoble and Lancaster 3 He experiments). Of these 
two procedures the first is more direct - vorticity formed 
during the quench is forced to migrate to the center of 
the container, through the existence of a subcritical rota- 
tion velocity, where it is detected. This permits loops of 
string of length larger than a known threshold to survive 
decay and results naturally in a higher density, i.e. the 
defect density is measured effectively very shortly after 
the transition takes place and need not be limited the the 
"infinitely long string". As a consequence, much smaller 
values of / and larger string densities are measured than 
in the 4 He experimental setting. 

Both remaining 3 He experiments end at a region of the 
phase diagram far from the transition - 3 He being very 
much colder than in the Helsinki experiment. A lower 
effective temperature results in the effective absence of 
damping mechanisms which in turn leads to the preser- 
vation of even small vortex loops. Dissipation mecha- 
nisms rely on the presence of quasiparticles. Thus, when 
the medium is very cold, energy dissipation will slow ex- 
ponentially and vortices can be stabilized by a coherent 
flow resulting from their motion through the supcrfluid. 
We expect therefore the string population in all 3 He ex- 
periments to be mostly in the form of relatively small 
loops. In Helsinki, the largest loops are stabilized by the 
slow rotation of the whole system, and their density can 
be extrapolated to the smaller loops, leading to the total 
consistent with the Grenoble and Lancaster experiments, 
where - one may guess - all of the loops survive for a long 
time in the absence of dissipation. 

The lifetime of these loops is thus expected to be much 
longer than that of thermal loops formed at a quench 
through the A-line in 4 He. As a result the long time 
decay of vorticity may also be very different in these two 
cases as the former corresponds to an ensemble of moving 
loops, at relative distances much larger than their typical 
radius, but the latter contains strings of all sizes, where 
the mean distance between strings is comparable to their 
length. 

This conjectured picture, supported in part by numer- 
ical studies, leads to the conclusion that both experi- 
mental settings in 3 He should lead to a positive result, 
compatible with a relatively small value of / relevant for 
all loops (we get / ~ 4, when n ~ l/(£/) 2 is used to fit 
early data in Fig. 6), whereas in 4 He the smallness of the 



signal at the time of first measurement makes the detec- 
tion more difficult and at present below the sensitivity 
threshold. 



DISCUSSION 

The mechanism we described early on in this paper is 
based on the analysis of the behavior of the order param- 
eter ip. The order parameter is clearly a phenomenologi- 
cal entity and the equations that govern its evolution are 
approximate and in many cases postulated rather than 
derived. On the other hand the underlying physics is usu- 
ally very specific. It may, for instance, involve atoms of 
some particular isotope such as 4 He. Thus, in principle, 
one could formulate an exact microscopic theory of par- 
ticular second order phase transformations. However, in 
all of the experimentally accessible cases discussed above 
such a fundamental theory is simply too complicated to 
lead to useful conclusions. The superfluid transition in 
4 He is a good example: Strong interactions in 4 He make 
it impossible to proceed rigorously all the way starting 
at the microscopic level. Analysis of related issues in 
the field theoretic context is also difficult [3^|. Recently 
however a new system has become experimentally ac- 
cessible: Atomic Bose-Einstein condensates (BEC's) un- 
dergo the second-order phase transition at much lower 
densities. Natural approximation schemes can be there- 
fore suggested, and the exact microscopic theory can be 
studied in greater detail than for the "old" superfluids. 

We shall not attempt to review the theoretical or ex- 
perimental situation in BEC's. Good reviews already 
exist (see eg. (36|). Our aim is simply to point out that 
questions concerning the formation of topological defects 
can be posed and analyzed within a much more funda- 
mental formalism, which is explicitly quantum. The ap- 
proximations start from the Schrodinger, equation and 
lead in a controlled manner to master equations for the 
density operator of the condensing system. Further ap- 
proximations result in a quantum kinetic theory. Prelim- 
inary analysis of these issues |37( allows one to recover 
the key scaling relations and the key predictions we have 
described in section II. Indeed, time dependent Landau- 
Ginzburg theory follows as an approximation to some of 
the terms which one obtains from the microscopic treat- 
ments. On the other hand, the microscopic theory con- 
tains additional terms, which alter predictions concer ning 
the formation of topological defects. Limited studies [37] 
indicate that the predicted densities of the vortex lines 
or of the winding numbers would be smaller than those 
based on the scalings of the order parameter (see section 
II of this paper). Moreover, corrections seem to be more 
significant as the ratio tq/tq decreases. 

The possibility of experimental studies of defect forma- 
tion in BEC quenches nevertheless exists and may lead 
to exciting insights into the problem. 
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Superconductors may be the other useful testing 
ground. Indeed, two experiments have been already re- 
ported [TH, with the claim of conflicting results, which 
seemed to depend on the geometry. Rapid cooling pro- 
duced no detectable signal in a high-temperature film, 
although it is far from clear whether any was expected. 
The original claim that the effect was ruled out at the 
"~ 10 3 level" was based on an overly optimistic predic- 
tion, which did not recognize that the total flux expected 
to arise in the experimentally studied geometry scales as 
n 1 / 4 , i.e. only with the fourth root of the total number 
of defects (rather than with the square root; see section 



4 of Ref. 2l| for discussion). 

The revised prediction is close to the claimed sensi- 
tivity of the experiment, and given the uncertainties in 
the critical exponents of the high-temperature supercon- 
ductor, as well as the possibility of imperfect trapping of 
the defects, etc. it is unfortunately impossible to extract 
useful constraints from the existing negative experiment. 

The experiment carried out by the same group, in the 
loop geometry has, on the other hand, yielded positive 
results. This experiment also operates near the edge of 
detectability. It detects the flux induced by a loop which 
is artificially broken into a large number N of super- 
conducting sections, which are then rapidly reconnected. 
The predicted flux should have a Gaussian distribution 
with a random direction and intensity corresponding to 
a ~ \f~N flux quanta [111, such signals appear to have 
been indeed found [18]. 

The available experimental results can be therefore de- 
scribed as confusing. In liquid crystals the results seem 
to be the perhaps least ambiguous, but they concern a 
(weakly) first order transition. In superfluids, 3 He is still 
the strongest case for the mechanism, especially since 
the results between all the experiments (carried out in 
quite different parameter regimes, and using very differ- 
ent techniques) are consistent. On the other hand the 
relevance of the Helsinki experiment for the cosmological 
scenario has been recently questioned by numerical ex- 
periments [38] which indicate that the vorticity generated 
in such settings may be induced by the flow imposed in 
the Helsinki experiment to facilitate the process of their 
detection. These simulations were carried out under a 
very idealized set of assumptions (which included very 
small fluctuations and axial symmetry) and conclusions 
of Ref. [38^ appear to be inconsistent with the experiment 
[3^ ]. but much more remains to be done to clarify this 
issue. Indeed, a fully 3-D study with large fluctuations 
is under way Further experimental and numerical 

studies to investigate the role of rotation in stabilizing 
vortex loops and to explore the implications for the A-B 
phase transition, etc., are nevertheless essential. 

The existing i He data are clearly disappointing, but 
not at odds with a more conservative theoretical esti- 
mate. Moreover as we have argued above there may be a 
way to reconcile estimates of vortex line density obtained 



from 3 He and 4 He experiments, even without any spe- 
cial appeal to the Ginzburg regime [32] or to quantum 
kinetic theory. 

Finally the first experimental reconnaissance into 
quenches in superconductors is preliminary in its nature 
and ambiguous in its results. 

In the meantime, numerical studies have confirmed and 
refined the basic indications of the theory of order param- 
eter dynamics. 

One may be therefore justified in the expectation of an 
exciting but uncertain future. A lot is at stake, includ- 
ing the understanding of the phase transition dynamics, 
nature of the order parameter and other collective ob- 
servables of quantum many body systems and perhaps 
even the relation between the quantum and the classical. 
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